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Abstract 

We propose an intensity-dependent pion-nucleon coupling Hamiltonian within a 
unitary multiparticle-production model of the Auerbach-Avin-Blankenbecler-Sugar 
(AABS) type in which the pion field is represented by the thermal-density ma¬ 
trix. Using this Hamiltonian, we explain the appearance of the negative-binomial 
(NB) distribution for pions and the well-known empirical relation, the so-called 
Wroblewski relation, in which the dispersion D of the pion-multiplicity distribution 
is linearly related to the average multiplicity < n > : D = A < n > +B, with the 
coefficient A related to the vacuum energy of the pion field and < 1 . The Hamil¬ 
tonian of our model is expressed linearly in terms of the generators of the SU{1, 1) 
group. We also find the generating function for the pion field. All higher-order 
thermal moments can be calculated from this function. At T = 0 it reduces to the 
generating function of the NB distribution. 
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1 Introduction 


During the last years a considerable amount of experimental information has been 
accumulated on multiplicity distributions of charged particles produced in pp and 
pp collisions in the centre- of-mass energy range from lOGeV to ISOOGeV . Mea¬ 
surements in the regime of several hundred GeV [1] have shown the violation of 
the Koba-Nielsen-Olesen (KNO) scaling [2], which was previously observed in the 
ISR c.m.energy range from 11 to 63 GeV [3]. The violation of the KNO scaling is 
characterized by an enhancement of high-multiplicity events leading to a broadening 
of the multiplicity distribution with energy. 

The shape of the multiplicity distribution may be described either by its 
C moments, Gq = (n'^)/(n)'^, or by its central moments (higher-order dispersions), 
Dq = ((n—g = 2, 3 .... The exact KNO scaling implies that all Gq moments 
are energy independent. Only at energies below lOOGeK do the G moments appear 
to be energy independent. It can also been shown [4] that the KNO scaling leads to 
a generalized Wroblewski relation [5] 

Dq = Aq{n) - Bq, (1) 

with the energy-independent coefficients Aq and Bq. The pp and pp inelastic data 
below lOOGeK also show the linear dependence of the dispersion on the average num¬ 
ber of charged particles, but with the coefficients Aq and Bq that are approximately 
equal within errors. 

The fact that the dispersion of the multiplicity distribution grows linearly with 
(n) implies that the elementary Poisson distribution resulting from the independent 
emission of particles is ruled out. 

The total multiplicity distribution of charged particles for a wide range of 

energies (22 — 900GeV) is found to be well described by a negative-binomial (NB) 
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distribution [1,6] that belongs to a large class of compound Poisson distributions [7]. 
The compound Poisson distribution is fully determined by its generating function 
G{z) of the form 

G{z) = ^ z^Pn, (2) 

where N is the average number of N independently produced clusters or clans (Pois¬ 
son distributed) [8]. Each of these subsequently decays according to the probability 
distribution corresponding to the generating function g{z). The NB distribution is 
obtained by choosing g{z) = ln{l — pz)/ln{l — p), with p = (1 — {n)/k)~^. It is a 
two-step process [8] with two free parameters: the average number of charged par¬ 
ticles (n) and the parameter k which affects the shape (width) of the distribution. 
The parameter k is also related to the dispersion H = ZI 2 by the relation 


D _.2 _ 1 ^ 

(n) k^ (n) ’ 

so that the observed broadening of the normalized multiplicity distribution with 
increasing energy implies a decrease of the parameter k with energy. The KNO 
scaling requires constant k. 


Although the NB distribution gives information on the structure of correlation 
functions in multiparticle production, the question still remains whether its clan- 
structure interpretation is simply a new parametrization of the data or has a deeper 
physical insight [9]. Measurements of multiplicity distributions in pp collisions at 
TeV energies [10] have recently shown that their shape is clearly different from that 
of the NB distribution. The distributions display the so-called medium-multiplicity 
’’shoulder”, with a shape qualitatively similar to that of the UA5 900GeV and UAl 
distributions [11]. A satisfactory explanation of this effect is still lacking [12]. 

In this paper we propose another approach to multiplicity distributions based 

on a unitary eikonal model with a pion-held thermal-density operator given in terms 
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of an effective intensity-dependent pion-nucleon coupling Hamiltonian. We assume 
that the system of produced hadronic matter (pious) is in thermal equilibrium at 
the temperature T immediately after the collision. 

The paper is organized as follows. In Sect.2 we explain the basic ideas of 
our unitary eikonal model with a pion-held thermal-density operator. A discussion 
of the Wroblewski relation is presented in Sect. 3. Finally, in Sect. 4 we draw 
conclusions and make remarks on the possible extension of the model to include 
two-pion correlations in the effective pion-nucleon Hamiltonian. 
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2 Description of the model 


At present accelerator energies the number of secondary particles (mostly pions) 
produced in hadron-hadron collisions is large enough, so that the statistical ap¬ 
proach to particle production becomes reasonable . Most of the properties of pions 
produced in high-energy hadron-hadron collisions can be expressed simply in terms 
of a pion-held density operator. We neglect difficulties associated with isospin and 
only consider the production of isoscalar ’’pions”. We expect that in high-energy 
collisions most of the pions are produced in the central region. In this region the 
energy - momentum conservation has a minor effect if the transverse momenta of 
the pions are limited by the dynamics and their rapidities restricted to the range 
\y \<Y, where Y = ln{s/m?‘) is the relative rapidity of the colliding particles. 

A. The AABS model 

A long time ago Auerbach, Avin, Blankenbecler, and Sugar [13] presented 
a class of models (AABS models) for which the scattering operator satished the 
exact s-channel unitarity at high energies. In these models, the incident hadrons 
propagate through the interaction region, without making significant changes in 
their longitudinal momenta (leading-particle effect). Only the part W = K^/s of 
the total c.m. energy ^/s in every concrete event is avaliable for particle production, 
where K is the inelasticity: 0 < A < 1 . 

In the AABS type of models, the scattering operator S is diagonal in the 
rapidity difference Y and in the relative impact parameter B of the two incident 
hadrons. The initial-state vector for the pion held is SiY^B) \ 0 >, where the 
vacuum state | 0 > for pions is in fact a state containing two incident hadrons. 

The n-pion production amplitude for n > 1 is given by 

*T„(T, B-h... K) = 2s{kr ...K\ S{Y, A) | 0). 
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We write the square of the n-pion production amplitude in the form 


I Tn{Y, B-h... kn) r= 4s2Tr{p(F, B) \ ki... kn){ki ... fcj}, (5) 

where the pion-density operator p(y, B) is dehned as 

p(y,B) = yr,B)|o>{o|St(y,B). (6) 

The square of the elastic scattering amplitude is then the matrix element of p{Y, B) 
between the states with no pions, i.e., (0 | p{Y, B) | 0). In terms of the pion-number 
operator 

^ H H ^ (7) 

k k 

the square of the ^'-matrix element when no pions are emitted can also be written 
in the form 

I (0 I ^(F, B) I 0) |2= Tr{p{Y, B) : :} = (8) 

Here :: indicates the operation of normal ordering and Q{Y,B) is the usual eikonal 
function (or the opacity function) of the geometrical model [14], The connection 
with the inelastic cross section and the exclusive cross section for production of n 
pions is then 

a,neiiY,B) = l-e-^^^'^\ ( 9 ) 

and for n > 1, it is 

aniY, B) = Tr{p(F, B) : — (10) 

In terms of a normalized pion-multiplicity distribution at each impact parameter, 
Pn(Y,B) = an(Y, B)/ainei(Y, B), the observed complete multiplicity distribution 
Pn{Y) is obtained by summing Pn{Y,B) over all impact parameters B with the 
weight function Q{Y, B) = ainei{Y, B)/(7inei{Y), i.e., 

Pn{Y) = J d^BQ{Y, B)Pn{Y, B). ( 11 ) 
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The first-order moment of PniX) is the average multiplicity 

(n) = Y.^Pn{Y) = J d^BQ{Y,B)n{Y,B). (12) 

The higher-order moments of Pn{Y) give information on the dynamical fluctuations 
from (n) and also on the multiparticle correlations. All these higher-order moments 
can be obtained from the pion-generating function 

G{z) = Y.z^Pn{Y) = I d‘^BQ{Y,B)G{Y,B;z), (13) 

by differentiation where 

G{Y,B-,z)=Tr{p{Y,B)z^} (14) 


is the pion-generating function in 5-space . Thus the normalized factorial moments 
Fg are 


F,= 


{n{n -l)...{n-q + l)) _ ,^,_^ d^G(l) 


(n)^ 


= {ny 


dz<i 


(15) 


and the normalized cumulant moments Kg are 


These moments are related to each other by the formula 

=!:('' 7 ( 17 ) 

For the Poisson distribution, all the normalized factorial moments are identically 
equal to 1 and all cumulants vanish for g > 1. 


We are concerned here mostly with the q = 2 moments, which are directly 
related to the dispersion D: 

F 2 = K 2 + I = (-rw)^ + 1 - ww- (18) 

(n) (n) 
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B. Thermal-density operator for the pion field 

The operator | 0)(0 | appearing in the definition of p{Y,B) represents the 
density operator p{vac) for the pion-field vacuum state. It can also be considered as 
describing the pion system in thermal equilibrium at the temperature T = 0. The 
density operator for a pion field in thermal equilibrium at the temperature T is then 

where 

Hq = ^ ct;fc(a|.afc + A), (20) 

k 

InZ = 

k k 

The quantity A '^k represents the lowest possible energy of the pion system at 
the temperature T = 0. The ” zero-point energy ” corresponds to A = A. If the 
pion energies Uk = \/k‘^ + are closely spaced, the summation over k is replaced 
by an integral: Xfc / cPk/2uJk- 

Note that p{vac) = Pt=o- The mean number of thermal (chaotic) pions is 

Ut = 


The transformed thermal-density 
Pt{Y,B) = 


E 


1 




( 21 ) 


E ^Tk- 


operator pt(Y,B) is now 


S{Y,B)ptS^Y,B) 




( 22 ) 


where 


H{Y,B)=S{Y,B)HoS^Y,B) 


( 23 ) 



is regarded as an effective Hamiltonian describing the pion system in the presence of 
two leading particles (nucleons). Taking into account an old observation of Golab- 
Meyer and Ruijgrok [15] that the Wroblewski relation can be satisfied for all energies 
if the square of the pion-nucleon coupling constant increases linearly with the mean 
number of pions (n), we propose the following form of the effective pion-nucleon 
intensity-dependent coupling Hamiltonian: 

H{Y,B) = ^[efc(F, B){Nu + A) + gu{Y, B){au ^+ 2A - 1 + h.c.)] (24) 

k 

= Y.Hk{Y,B), 

k 

where e|(H, B) = ul + AgliY, B). The interaction part of the Hamiltonian for the 
k mode is no longer linear in the pion-field variables and represents an intensity- 
dependent coupling [16]. It is also easy to see that the operators 

Ko{k) = Nk + X, (25) 

K-{k) = ak\JNk -f 2A — 1, 

K^[k) = Nk 2A — 1 Ofc 

form the standard Holstein-Primakoff [17] realizations of the sm( 1, 1) Lie algebra, 
the Casimir operator of which is 

Ck = Kl{k) - ]^[K+{k)K_{k) + K_{k)K+{k)] = A(A - 1)4. (26) 

The Hamiltonian Hk{Y,B) = is thus a linear combination of the generators of 
the SU (1,1) group: 

Hk = eMk) + gk[K+{k) + K_{k)]. (27) 

The corresponding S-matrix which diagonalizes the Hamiltonian H{Y, B) is 

Ars) = n&(rs), ( 28 ) 

k 
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where 

Sk{Y,B) = exp{-dk{Y,B)[K^{k) - K_{k)]}, ( 29 ) 

with 

th4(y,^) = (30) 

ek{Y,B) 

Since the dependence on the variables Y, B is contained only in the hyperbolic angle 
OkiY, B), we shall from now on assume this dependence whenever we write 6^.. 

It is easy to see that the initial-state vector for the pion held, S(Y,B) \ 0), 
factorizes 

•S(rs) |0) = n(St(VB) |0i.)), (31) 

k 

with 

S,(K.B)|0») = (1 I nO (32) 

nfc!r(2A) 

= \eu). 

where | Uk) = (nfc!)“^'^^(a[)"''= | 0^). In the same way we hnd that the pion thermal- 
density operator pt(Y,B) is also factorized as 

PT{Y,B) = l[pT{ek), (33) 

k 

with 

Pt(4) = ^ E I rik, 9k){nk, Ok |, (34) 

^k Uk 

where | nk,9k) = Sk(Y, B) \ Uk). Note that | nk,9k) form a complete orthonormal 
set of eigenvectors of the /c-mode Hamiltonian Hk, i.e., 

Hk \ nk,9k) = Uk{nk + \) \ nk,9k). (35) 
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3 Pion-generating function and its moments 


The average multiplicity nT{Y,B), the dispersion and all higher-order 

moments 

nUY, B) = Tr{pT{Y, B)^}, g = 1, 2,... (36) 

at the temperature T in i? space can be obtained from the pion-generating function 


GT{Y,B-z) = X{GT{du,z) 

k 


by differentiation, where 


Grief,z) ^Tr[pT(ek)z”‘‘}. 


(37) 


(38) 


After performing a certain amount of straightforward algebraic manipulations, 
we arrive [18] at the following expression for the pion-generating function GT{ 0 k', z): 

GriOk; z) = Go{e,; ;^)(1 - R^\l + y, + (39) 

where 

Rk = \/l - 2xkyk + yl, (40) 

_ z + { 1 - zfsh'^{9k)di^{9k) 
z - {I - zYsh‘^{dk)ch‘^{9ky 

^ 0 PP, z - {I - z)sh\dk) 
i + ii-zYshyek)' 

and Go(^fc; z) denotes the pion-generating function at the temperature T = 0: 

Go^Ok, ;2) = [1 -f (1 - z)sh^{ek)]~‘^^. (41) 

We observe that Gq is exactly the generating function of the NB distribution with a 
constant shape parameter 2A, and the average number of k-mode pions is equal to 

niOk) = 2 \sh^{ 6 k). 
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The vacuum value of the /c-mode thermal-density operator pxiOk) is used to obtain 
the /c-mode thermal eikonal function 


(Ofc I pr(4) I Ofc) = 




The total eikonal function is flriX, B) = Y.k ^riOk)- 

For the /c-mode pion held in B space at the temperature T, we hnd the following 
average number and the dispersion: 


nriOk) = n{ek) riTk + ^n(6'fc)nrfc, 

2A — 3 4 

d%{9k) = c0p^ + (f‘{9k)[^-\ -—h-Tfe + 


where 


= ^Tk + ^Tk-, 

= ^fi^{9k)+n{9k) 


Two limiting cases are of interest namely, T = 0 and T —cx). For the T = 0 


case, we hnd 


d\9k) _ 1 1 

n^{9k) 2A h(6'fc)’ 


as it is to be expected from the NB distribntion ( Eq. 3). However, our interpretation 
of this result is quite diherent. In onr case, the parameter A is connected with the 
vacnnm energy of the pion held in Eq. 20, and has nothing to do with either the 
nnmber of pion sources or the number of clans. It is a real and positive constant 
labeling the positive discrete class of a unitary irredncible representation of SU (1,1), 
which is a dynamical symmetry gronp of onr system Hamiltonian. It is important 
to observe that pions in the k mode are distributed according to the NB distribution 



with a constant shape parameter 2A. The Wroblewski relation 


d(et) = An(»,) + B 


( 47 ) 


is obtained with energy-independent coefficients A = (2A) and B = (A/2)^/^. If 
A > 1/2, we have A < 1. 


The contribution from all the k modes gives 


(f{Y,B) 

n^{Y,B) 




1 

h(F, B) ’ 


(48) 


where p{0k) = n{6k)/n{Y, B). In this case, the coefficient A in the Wroblewski 
relation becomes energy dependent and is of the form 


k 

Since J2kP{^k) = 1 and all p{9k) are positive functions of 9k, 
always smaller than one. Therefore, A{Y,B) < 1 if A > 1/2. 


( 49 ) 

the sum J2kP‘^{9k) is 


Finally, the summation over all impact parameters gives 

(/l)^ = /</^BQ(y,B)[(.4^(y.B) + l)(d|^)2-l| + i. (50) 

This expression, when combined with our preceding analysis suggests that the coef- 
hcient A in the Wroblewski relation should be energy dependent and smaller than 
one. 


For the temperature T going to inhnity we obtain 


dU9k 


n\k{9k 


T^OO 


= 2-(l + 


n{9, 


k )\-2 


A 


= l + th^{29k). 


(51) 


This result shows that at very high temperature the distribution of pions will become 
chaotic if 9k is very small. This will happen when ujk 3> gk{Y,B). 


13 



4 Conclusions 


In this paper we have proposed an intensity-dependent pion-nncleon coupling Hamil¬ 
tonian with SU{1, 1) dynamical symmetry, within a multiparticle-production model 
of the AABS type in which the k mode pion held is represented by the thermal- 
density operator. We have shown that this Hamiltonian explains in a natural way 
the appearance of the NB multiplicity distribution for pions in impact-parameter 
space.The shape parameter of the NB distribution is related to the vacuum energy 
of the pion held at the temperature T = 0. 

The Wroblewski relation is obtained with the coefficient A that is energy de¬ 
pendent and smaller than one if the vacuum energy of the pion held is larger than 
’’zero-point energy” corresponding to A = 1/2. 

For T 7 ^ 0, we have found a pion-generating function that may be used for 
obtaining all higher-order moments of the pion held. 

In our model, the k modes of the pion held are statistically independent and are 
described by the factorized thermal-density operator. Correlations between diherent 
k modes are absent and, at this stage, our model cannot describe the emission of 
resonances. However, this can be remedied by adding a mode-mode interacting part 
to the total Hamiltonian H(Y,B), e.g., J2k,k'[-^k,k'(X, + h.c.] [19]. We hope 

to treat this case elsewhere. 
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